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摘 要 ;水 泥土 材料 在 剪 切 破坏 过 程 中 常会 产生 剪 胀 现象 ,但 目前 基于 Mohr-Coulomb 强度 理论 的 
恒定 剪 胀 角 模 型 不 能 正确 表达 水 泥土 材料 在 变形 破坏 过 程 中 的 非 线 性 体积 变化 行为 。 因 此 ,根据 
不 同 水 泥 掺 入 比 和 围 压条 件 下 的 水 泥土 三 轴 压 缩 试验 ,结合 塑性 力学 理论 ,采用 非 线 性 拟 合 方法 建 
立 了 能 同时 考虑 单 轴 抗 压强 度 、 平 均 应 力 和 等 效 塑 性 应 变 影响 的 剪 胀 角 模 型 ,并 运用 有 限 元 技术 进 
行 了 数值 模拟 验证 。 结 果 表 明 ,水 泥土 三 轴 试 验 的 应 力 - 应 变 曲 线 呈 软化 型 ,水 泥土 发 生 破 坏 后 会 
伴随 着 明显 的 体积 膨胀 。 水 泥土 的 剪 胀 行为 通常 用 剪 胀 角 来 量化 ,考虑 平均 应 力 和 等 效 塑性 应 变 
影响 的 前 胀 角 模型 能 够 合理 表达 水 泥土 在 峰 后 的 非 线性 体积 变化 行为 ,模拟 结果 与 试验 结果 具有 
一 较 好 的 一 致 性 。 而 将 剪 胀 角 取 为 恒定 值 会 出 现 高 估 或 低估 剪 胀 效应 的 情况 ,不 能 正确 反映 水 泥土 
峰 后 的 非 线性 体积 变化 规律 。 
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Triaxial test and numerical simulation verification of dilatancy 
characteristics of cemented soils 
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Abstract :The dilatancy usually occurs in cemented soils materials when it is at failure by shear, but the 
constant dilatancy angle model based on Mohr-Coulomb strength theory at present cannot correctly express 
the nonlinear volume change behavior of cemented soils materials during deformation and failure. There- 
fore , according to the triaxial test of cemented soils under different cement ratios and confining pressure 
conditions , combined with the theory of plastic mechanics , the model of dilatancy angle affected by the uni- 
axial compressive strength , equivalent pressure stress and equivalent plastic strain simultaneously is estab- 
lished with a nonlinear fitting method , and then , the model of the dilatancy angle is verified by numerical 
simulation using finite element technology. The results show that the stress-strain curves of the cemented 
soils triaxial test is a softening type, and an obvious volume expansion will be generated when the cemented 


soils is damaged. The dilatancy behavior of cemented soils is usually quantified by the dilatancy angle. The 
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dilatancy angle model under the influence of equivalent pressure stress and equivalent plastic strain can 


reasonably express the nonlinear volume change behavior of cemented soils after the peak stress , and the 


simulation results are in good agreement with the test results. However, if the dilatancy angle is taken as a 


constant value, the dilatancy effect will be overestimated or underestimated , and it cannot correctly reflect 


the law of nonlinear volume change of cemented soils after the peak stress. 


Key words :cemented soil ;triaxial test ; volume expansion; dilatancy angle model 


在 建设 工程 中 经 常会 遇见 软 土地 质 条 件 ,在 这 
些 软 土地 区 进行 如 修筑 公路 建筑 基 坑 开 挖 河 岸 加 
固 等 工程 时 ,通常 需要 对 软 土 进行 加 固 处 理 。 水 泥 
土 桩 是 加 固 软 土 的 常用 技术 , 常 采用 深层 搅拌 或 高 
压 喷射 注 浆 施 工 将 水 泥浆 液 与 原 土 混合 ,以 实现 对 
软 土 的 加 固 。 它 可 以 最 大 限度 地 利用 地 基 土 ,使 原 
有 地 基 的 强度 和 刚度 得 到 较 大 提高 , 且 具 有 造价 低 
厥 工期 较 短 等 优点 ,因而 在 公路 .铁路 水利, 市 政 
以 网 建筑 等 工程 建设 中 得 到 了 广泛 的 应 用 "” 。 
水 泥土 相对 于 原 土 ,在 强度 特性 、 变 形 特性 、 渗 
透 竺 性 和 动力 特性 等 方面 都 有 着 较 大 的 改善 ,针对 
不 同 的 工程 应 用 目的 ,水 泥土 被 关注 的 工程 特性 也 
不 同 。 水 泥土 的 强度 与 变形 特性 是 被 关注 较 多 的 地 
7 六 两 个 性 质 影响 着 被 加 固 对 象 的 工程 性 能 ,如 地 
基 的 承载 力 和 稳定 性 5 , 基 坑 的 土 压 力 和 变形 , 
软 韦 边 坡 的 稳定 性 和 变形 等 ”。 然 而 , 随 着 水 泥土 
材料 在 工程 中 的 广泛 应 用 ,其 工作 条 件 及 所 处 的 应 
为 器 态 日 益 复杂 ,常规 的 设计 计算 已 无 法 满足 工程 
建设 的 需要 ,全 面 认识 水 泥土 的 强度 特性 和 本 构 模 
型 6 夫 而 开展 数值 模拟 是 当前 对 水 泥土 应 用 的 研究 
超 棚 '””) 。 众 多 研究 者 发 现 ,水 泥土 应 力 -应 变 曲线 
表 钢 出 超 固 结 土 的 特性 ,强度 经 过 峰值 后 有 软化 特 
征 , 剪 切 变形 过 程 有 剪 胀 现象 。LO 等 "通过 对 粉 
土 的 三 轴 试 验 发 现 ,经 水 泥 混合 后 的 粉 土 在 破坏 时 
相 比 原 土 产生 了 较 大 的 前 胀 ,前 胀 量 可 以 用 剪 胀 角 
表征 , 当 水 泥 挫 量 增加 ,水 泥土 变 得 越 硬 脆 , 剪 胀 角 
也 越 大 , 而 围 压 则 产生 相反 的 关系 。SARIOSSEIRI 
等 ' 从 水 泥土 的 不 排水 三 轴 试 验 结果 得 出 ,经 水 泥 
处 理 后 土 的 抗 剪 强度 显著 提高 ,但 破坏 类 型 也 发 生 
了 较 大 的 改变 , 随 着 水 泥 挨 量 的 增加 , 土 的 破坏 由 延 
性 破坏 向 脆性 破坏 过 渡 , 体 现在 水 泥土 破坏 时 的 应 
变 由 4% 降 为 略 大 于 1% ,同时 伴随 着 试 样 体积 的 前 
切 膨 胀 。 宋 新 江 等 '” 通过 轴 对 称 条 件 下 水 泥土 强 
度 特性 试验 研究 ,发现 剪 切 初始 阶段 ,水 泥土 表现 出 
一 定 的 剪 缩 。 随 着 剪 切 强度 的 进一步 加 大 ,水 泥土 
各 连接 结构 产生 破坏 ,颗粒 之 间 发 生 相互 错 动 ,体积 


逐渐 增加 ,表现 出 一 定 的 剪 胀 。CHADAKPOUR 
等 ”研究 了 黄土 和 黏土 的 三 轴 剪 切 力学 行为 ,发 现 
水 泥 能 显著 改善 条 土 的 力学 性 能 ,使 试 样 由 可 塑性 
向 脆性 发 展 。 水 泥 的 掺 和 使 试 样 产生 了 剪 胀 行为 ， 
并 受到 水 泥 掺 量 和 围 压 的 影响 。 

在 连续 介质 理论 中 ,最 广泛 用 来 衡量 岩 土 材料 
扩容 和 控制 材料 体积 变化 的 参数 是 剪 胀 角 "" 。 然 
而 ,关于 水 泥土 在 宕 土工 程 中 的 应 用 计算 , 剪 胀 角 经 
常 被 忽视 ,即使 考虑 ,通常 采用 假设 剪 胀 角 为 0 或 其 
他 恒定 值 的 方式 进行 处 理 “ 。 剪 胀 角 的 简化 假设 ， 
无 法 真实 反映 水 泥土 的 变形 破坏 过 程 ,降低 了 对 工 
程 结构 性 能 判定 的 准确 性 。 材 料 的 剪 胀 变化 行为 受 
到 围 压 及 塑性 应 变 发 展 的 影响 " ,水 泥土 的 剪 胀 行 
为 在 已 有 研究 中 已 经 被 证 明 受 多 因素 的 影响 ,前 胀 
角 的 变化 是 一 个 非 线 性 的 过 程 ,但 正确 识别 水 泥土 
的 影响 因素 .并 建立 起 剪 胀 角 的 演化 模型 这 方面 的 
工作 还 鲜 有 研究。 因此 ,本 研究 通过 三 轴 压 缩 试验 ， 
研究 水 泥土 的 剪 胀 特性 ,并 建立 起 剪 胀 角 的 演化 模 
型 ,最 后 通过 数值 模拟 对 所 建立 的 水 泥土 剪 胀 模型 
进行 验证 ,以 期 为 相关 工程 提供 参考 。 


c— 


l 试验 方案 


1.1 试验 材料 


本 次 试验 用 土 取 自 广西 梧州 长 洲 水 利 枢纽 三 线 
四 线 船闸 工程 右岸 边 坡 的 一 、 二 级 阶地 冲积 层 。 所 
取 试验 用 士 为 灰白 色 粉 质 黏 士 , 按 《 土 工 试验 方法 标 
ME) (GB/T 50123 - 2019) 对 试验 土 样 的 物理 力学 性 
质 指标 进行 测定 ,结果 见 表 1。 采 用 密度 计 法 对 试 
验 土 样 进行 颗粒 分 析 , 结 果 见 图 1。 由 图 可 见 土 样 
颗粒 级 配 连续 。 试 验 所 用 水 泥 为 华润 P. 0 42.5 级 
普通 硅 酸 盐水 泥 , 水 泥水 经 过 与 土 混合 并 固化 后 形 
成 水 泥土 。 


1.2 试验 内 容 及 方法 
本 研究 主要 从 水 泥土 的 水 泥 掺 量 、 塑 性 变形 和 
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HIE 3 个 影响 因素 作用 下 去 测量 水 泥土 剪 胀 角 的 变 
化 规律 。 为 了 较 清 楚 地 区 分 水 泥土 中 水 泥 成 分 和 原 
土 成 分 的 关系 ,以 反映 水 泥 用 量 大 小 ,定义 水 泥 摊 入 
EK: KEBA Ha.) = (水 泥 质 量 / 干 土质 量 ) x 
100% 。 

为 了 研究 水 泥 摊 量 对 水 泥土 前 胀 性 的 影响 , 设 
置 了 4 个 不 同 的 水 泥 摊 入 比 (10%、20% 、30% 、 
4096 ) ,这 些 水 泥 挫 量 覆 盖 了 水 泥土 工程 应 用 中 可 能 


围 是 100 ~ 500 kPa ,这 种 压力 状态 能 够 代表 水 泥土 
结构 在 实际 工程 中 可 能 受到 的 应 力 ! … 。 由 于 所 研 
究 问题 针对 地 基 处 理 , 水 泥土 基本 处 于 饱和 状态 ,所 
以 不 考虑 含水 量 的 影响 ,在 制作 试 样 时 加 入 相同 的 
用 水 量 ,并 用 浸水 养护 方式 对 试 样 进行 人 饱和。 在 选 
择 三 轴 试 验 的 类 型 时 ,考虑 到 水 泥土 刚度 大 、 变 形 
小 、 含 水 量 低 、 孔 际 水 压力 相对 于 水 泥土 的 抗 压强 度 
来 说 微小 ,所 以 整个 试验 过 程 忽 略 孔 际 水 的 作用 , 采 


使 用 到 的 水 泥 用 量 。 三 轴 压 缩 试验 中 选取 的 压力 范 。 用 的 是 固 结 排水 前 切 试验 ( CD)。 
表 1 试验 土 样 的 物理 力学 性 质 指 标 
Tab.1 Material properties of soil samples for test 
含水 量 P mm 饱和 度  0RBR 塑 限 塑性 Werk RRI 内 摩擦 角 ”压缩 系数 
EEZ P w/ á E E Si/ pw wrio/ Wp/ HA HA c/ w/ A 
p See s % > k % h h, kPa — (?) MPa! 

gE 22.5 1.99 0.690 89.10 2.740 34.29 19.10 15.19 0.22 20 17 0.29 
PR 

C 轴 仪 , 试 样 在 规定 围 压条 件 下 固 结 完成 后 再 进行 加 
e 载 ,控制 加 载 速率 为 0.6 mm/min。 试 验 数据 由 电脑 
d 自动 采集 ,采样 步 长 0.2 mm. 

= 2 试验 结果 分 析 

2 2.1 水 泥土 变形 破坏 及 剪 胀 特征 

= 水 泥土 三 轴 试验 得 到 的 应 力 -应 变 、 体 积 应 变 - 
"S 粒 径 /mm 轴 向 应 变 关系 曲线 如 图 2 所 示 。 

c 图 1 土 的 粒 径 级 配 累积 曲线 从 图 中 可 以 看 出 ,水 泥土 的 应 力 -应 变 关系 总 体 
C Fig. 1 Grain size distribution of the soil 呈 应 变 软 化 型 ,从 水 泥土 的 变形 破坏 过 程 可 见 ,水泥 
思 制 样 时 ,将 从 工程 现场 取得 的 土 料 经 烘 干 . 碾 。” 土 材料 具有 类 似 混凝土 这 种 脆性 材料 的 变形 破坏 特 


WE x2 mm 得 后 备用 。 拌 和 水 泥土 用 水 量 全 部 统一 
为 26% , 即 :w = 水 的 质量 /( 水 泥 质 量 + 干 土质 
Æ) =26% 。 为 保证 混合 料 的 均匀 , 先 将 干 土 和 水 泥 
粉 在 搅拌 倪 内 均匀 拌和 ,然后 再 通过 喷 水 器 往 倪 内 
喷洒 水 花 , 并 同时 搅拌 水 泥土 。 达 到 搅拌 用 水 量 要 
求 后 ,将 搅拌 均匀 的 水 泥土 装 入 直径 为 39.1 mm \ 高 
度 为 80 mm 圆柱 体 三 为 模 中 制 样 。 由 于 水 泥土 呈 软 
塑 状 ,只 需 用 击 锤 分 5 层 捣 实 即 可 。 经 测试 ,所 制 水 
泥土 试 样 密 度 均 匀 ,平均 密度 为 2.028 4 g/cm , 因 拌 
土 用 水 量 与 击 实 功 相同 ,确保 了 试 样 具有 相近 的 密 
实 度 。 为 了 模拟 现场 实际 情况 ,使 试 样 接近 饱和 状 
态 , 试 样 在 装置 内 项 置 48 h 终 凝 后 脱 模 ,然后 放 入 
水 箱 中 浸水 养护 28 d 后 开始 三 轴 压 缩 试验 。 三 轴 
试验 使 用 南京 土壤 仪器 厂 生 产 的 TSZ 系列 全 自动 三 
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征 。 从 图 中 水 泥土 的 体积 应 变 - 轴 向 应 变 曲 线 还 可 
看 出 ,水 泥土 试 样 发 生 破坏 后 产生 了 较 大 的 体积 膨 
胀 ,这 个 现象 同样 也 被 其 他 学 者 所 证 实 '"。 综 合 
来 看 ,可 以 把 水 泥土 的 应 力 - 应 变 曲线 和 体积 应 变 - 
轴 向 应 变 曲线 概括 为 如 图 3 所 示 的 4 个 阶段 。 

第 1 阶段 为 弹性 变形 阶段 (04 阶段 ) ,应 力 -应 
变 曲 线 近似 为 一 条 直线 ,主要 是 水 化 凝固 的 水 泥土 
颗粒 受 力 产 生 弹 性 变形 。 曲 线 初 始 段 为 下 凸 的 弯曲 
段 ,主要 是 垂直 于 压 应 力 方 向 的 微 裂 纹 `. 孔 洞 以 及 试 
样 和 仪器 间 的 接触 间 陀 受 压 闭合 而 产生 的 ,通常 忽 
略 不 记 。 直 线段 的 峰值 点 A 称 为 弹性 极限 ,对 应 的 
应 变 在 0.01 左右 ,弹性 阶段 的 极限 荷载 与 极限 抗 压 
强度 的 比值 大 致 在 0.5 ~0.8 之 间 变 化 。 在 这 一 阶 
段 中 体积 变形 处 于 可 以 恢复 的 弹性 压缩 状态 。 
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图 2 水 泥土 应 力 -应 变 、 体 积 应 变 - 轴 向 应 变 关系 曲线 
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Fig.2  Stress-strain curves and volumetric-axial strain curves of cemented soils 


第 2 阶段 为 塑性 硬化 阶段 (48 阶段 ) ,水 泥土 的 
塑性 变形 开始 增长 并 产生 应 变 硬化 现象 。 该 阶段 
中 , 试 样 出 现 裂纹 扩展 ,应 力 进 入 缓慢 增长 阶段 直到 
峰值 点 C, 该 点 称 为 水 泥土 的 极限 强度 。 极 限 强 度 
受 水 泥 挫 入 比 和 围 压 的 影响 ,表现 为 水 泥 掺 入 比 和 
围 压 的 增加 ,水 泥土 的 极限 强度 明显 增长 。 如 图 4 
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所 示 , 水 泥 摊 量 和 围 压 也 对 破坏 应 变 有 着 显著 影响 。 
随 着 水 泥 掺 入 比 的 增加 ,达到 峰值 应 力 时 的 破坏 应 
变 有 所 减 小 。 围 压 越 大 ,水 泥土 塑性 硬化 过 程 越 长 ， 
软化 效应 减弱 ,破坏 应 变 得 到 提高 。 该 阶段 水 泥土 
试 样 的 体积 变形 ,开始 由 压缩 向 剪 胀 转化 。 随 着 轴 
向 应 变 的 增加 ,体积 膨胀 速度 加 快 ,也 意味 着 裂纹 进 
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一 步 扩 大 , 当 应 力 至 峰值 点 8 时 ,总 的 体积 应 变 大 致 
位 于 0 的 左右 位 置 。 


o 
Y 


3 水 泥土 受 压 典型 应 力 -应 变 曲线 和 
轴 向 应 变 - 体 积 应 变 曲线 


Fig.3 Typical compression stress-strain curves and 


volumetric-axial strain curves of cemented soils 
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ig.4 Variation of failure strain versus confine stress 
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and ration for cemented soils 

第 3 阶段 为 峰值 后 的 软化 阶段 (BC 阶段 ) ,也 是 
水 泥土 主要 的 破坏 阶段 。 加 载 进 入 该 阶段 后 , 试 样 
内 部 的 裂纹 向 外 扩展 ,逐渐 发 展 为 纵向 裂纹 并 在 斜 
方向 连接 贯通 ,形成 较为 明显 的 剪 切 面 。 在 这 个 阶 
段 ,应 力 下 降 较 快 ,应 力 -应 变 曲 线 向 下 弯曲 出 现 拐 
点 C。 随 着 变形 的 增加 ,应 力 -应 变 曲线 凸 向 水 平 轴 
发 展 并 趋 于 稳定 。 在 水 泥土 贯通 的 裂 际 中 ,水 泥土 
破碎 成 较 大 的 颗粒 充填 于 其 中 , 剪 胀 效应 在 该 阶段 
继续 发 展 , 但 剪 胀 的 速率 逐渐 变 缓 。 

经 过 收敛 点 C 以 后 进入 第 4 阶段 , 即 残余 强度 
阶段 。 该 阶段 试 样 的 应 变 继续 增 大 , 试 样 中 间 形 成 
一 条 贯穿 斜 截面 的 主 裂 颖 ,也 有 些 形成 交叉 的 两 条 
斜 裂缝 , 受 力 均匀 的 试 样 还 会 形成 锥 形 压 裂 , 试 样 的 
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破坏 形状 如 图 5 所 示 。 残 余 强 度 由 破坏 面 上 的 摩 阻 
力 和 残余 的 粘 结 力 提供 ,残余 应 力 o, 与 峰值 应 力 or, 
的 比值 总 体 上 随 着 围 压 的 增 大 而 增 大 (图 6) ,表现 
出 向 理想 塑性 变形 方向 发 展 的 趋势 。 在 体积 变化 方 
面 ,虽然 水 泥土 进入 残余 强度 阶段 ,但 仍然 会 产生 剪 
胀 现象 , 剪 胀 速率 随 着 水 泥土 试 样 应 变 的 增加 呈 非 
线性 降低 趋势 。 


图 5 水 泥土 试 样 变形 破坏 图 
Fig.5 Failure modes of cemented soils specimens 
100 


y-0.069x427.278 
R=0.524 8 


(0./0,)/% 


Gy kPa 
6 水 泥土 0,/ 0-05 的 变化 趋势 


Fig.6 Variation of ou/o versus aa for cemented soils 


2.2 水 泥土 剪 胀 效 应 的 表征 及 演化 规律 


剪 胀 性 可 以 定义 为 材料 中 单元 剪 切 变形 引起 的 
体积 变化 。 用 剪 胀 角 来 描述 材料 剪 胀 行为 是 一 种 比 
较 合理 的 选择 ,因为 它 代 表 塑 性 体积 应 变 与 塑 
性 剪 切 应 变 的 比值 。 以 简单 的 直 剪 试验 为 例 , 剪 
效果 可 描述 为 图 7 所 示 形 状 , 则 剪 胀 角 消 为 
Ah _ dei 
AT dy' (1) 

在 Mohr-Coulomb 55 HE 3 ie rp , IU (1) 4E X. HJ B 
胀 角 可 在 强度 准则 的 破坏 包 线 上 给 出 更 直观 的 摘 
述 。 如 图 8 所 示 ,根据 Drucker 公设 ,材料 届 服 后 , 塑 
性 势 方向 垂直 于 破坏 包 线 ,塑性 剪 切 应 变 增 量 与 总 
的 塑性 应 变 增 量 成 y 的 夹 角 , 当 y= 时 , 则 认为 材 
料 服从 相关 联 的 流动 法 则 , 当 y 关 9 时 , 则 认为 材料 
服从 非 关 联 的 流动 法 则 ” 。 而 在 三 轴 试 验 中 , 剪 切 
应 变 无 法 直接 测 出 ,因此 VERMEER 45! 先 建立 塑 


tany = 
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性 势 函数 形式 ,通过 测量 与 塑性 势 函 数 相关 的 塑性 
应 变 推导 出 前 胀 角 的 表达 式 , 即 


de 
siny LAEN (2) 
E, 


图 7 AEE KR 


Fig.7  Dilatancy effect of geotechnical materials 


图 8  Mohr-Coulomb 强度 准则 塑性 势 方向 
Fig.8 Plastic potential direction of Mohr-Coulomb 


strength criterion 


加 由 式 (2) 可 知 ,计算 出 剪 胀 角 的 关键 在 于 求 出 


JEN de? 。 借 鉴 ALEJANO 等 ”所 给 出 的 方法 ,如 


XB vk NEAR 2 


图 9 剪 胀 角 模型 构建 图 解 
Fig.9 Schematic of constructing dilatancy angle model 
假设 加 载 、 印 载 过 程 弹性 模 量 E 为 一 常数 , 则 从 
图 中 可 知 
p 


e =E -el (3) 


a =g e, (4) 


因为 
€; = (0, -04,)/E (5) 
£; =81 +283 - & - 2vei (6) 
从 而 联 立 式 (3) ~ (6) 可 得 
ei -6e,- (0, -0,)/E (7) 
e ze,-(1-2v) (o, -0,)/E (8) 


其 中 弹性 模 量 & 可 通过 应 力 -应 变 曲线 的 弹性 阶段 

计算 获得 ,而 泊 松 比 v 也 可 由 弹性 阶段 的 压缩 体积 

应 变 与 平均 应 力 关 系 求 得 , 即 

1 Ee, 

2 2(o -on) (9) 
在 图 9 中 ,虚线 反映 的 是 a? - 6l 的 关系 , 按 式 

(2) 计算 剪 胀 角 , 则 需要 知道 a? - el 关系 曲线 上 每 

一 点 的 变量 增 量 。 实 验 得 到 的 et - et 关系 曲线 是 由 

离散 的 点 构成 ,不 能 直接 求 出 微分 增 量 ,因此 本 研究 

采用 中 心 差分 方法 计算 剪 胀 角 。 如 对 于 曲线 上 点 

的 剪 胀 角 ,相应 的 计算 公式 为 


o del 
amp; a -2del 


v 


e iui ER 

“(a ga) A a 

由 式 (10) 可 知 ,前 胀 角 受 材料 塑性 变形 的 影 

响 , 在 水 泥土 的 变形 破坏 过 程 中 会 发 生变 化 ,而 引起 

这 种 变化 的 影响 因素 则 是 一 个 复杂 的 问题 。 在 塑性 

理论 中 ,通常 将 影响 力学 参数 的 变量 定义 为 塑性 功 、 

等 效 塑 性 应 变 或 塑性 体积 应 变 等 。 因 此 ,本 研究 统 

一 定义 剪 胀 角 为 等 效 塑 性 应 变 er 的 函数 ,通过 建立 

Y- e! 的 关系 曲线 ,反映 前 胀 角 在 水 泥土 剪 切 破坏 过 

程 中 所 呈现 的 变化 规律 。 在 三 轴 试 验 条 件 下 ,等 效 
塑性 应 变 e? 计算 式 为 

P 42 


(10) 


aan et ce «Cet e! e Cet cet) 
=el opel (n) 


在 三 轴 试 验 的 应 力 -应 变 关系 曲线 中 ,计算 出 每 
一 个 数据 点 对 应 的 剪 胀 角 和 等 效 塑性 应 变 , 从 而 可 
建立 三 轴 试 验 的 y-e? 关系 散 点 图 ,如 图 10 所 示 。 
由 图 可 知 , 剪 胀 角 随 着 等 效 塑 性 应 变 的 增加 , 先 陡 增 
至 最 大 值 后 呈 非 线性 降低 , 当 等 效 塑 性 应 变 达 到 
0.1 时 ,部 分 试 样 的 剪 胀 角 可 降低 至 最 大 剪 胀 角 的 
1Z10。 初 始 陡 增 段 的 等 效 塑 性 应 变 微小 ,大 约 在 
0. 002 ~ 0. 005 之 间 , 相 对 于 整体 的 变形 量 来 说 , BE 
增 段 几 乎 可 以 忽略 不 记 。 在 所 选择 的 围 压 范围 内 ， 
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围 压 的 增加 会 降低 剪 胀 角 ,但 影响 没有 等 效 塑 性 应 
变 大 。 另 外 ,水 泥 挫 和 人 比 对 剪 胀 角 也 有 一 定 的 影响 ， 


但 总 体 影响 不 大 。 试 样 的 加 载 过 程 也 是 塑性 应 变 和 
30 
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图 10 水 泥土 三 轴 试 验 水 - e 关系 散 点 图 


应 力 状态 不 断 改变 的 过 程 ,同时 还 受 水 泥 掺 入 比 的 
影响 ,从 而 决定 了 剪 胀 角 与 等 效 塑 性 应 变 应 力 状态 
以 及 水 泥 挫 入 比 之 间 存 在 着 客观 上 的 联系 。 
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pa Fig. 10 Scatter diagram of y - er relationships of cemented soils using triaxial tests 


2 x 水 泥土 剪 胀 角 演 化 模型 构建 


因 剪 胀 角 的 变化 受到 等 效 塑 性 应 变 、 围 压 以 及 
AGRIS A tI E 所 以 , 剪 胀 角 模 型 可 以 表达 为 3 
个 饮 响 因素 的 函数 , 即 y = f(g?,g,,a.) 。 三 维 应 力 
状态 下 ,最 小 主 应 力 无 法 单方 面 反映 应 力 状 态 对 剪 
胀 角 的 作用 ,考虑 更 一 般 情况 下 ,用 平均 应 力 ww fX 


虑 ,最 终 剪 胀 角 模 型 方程 表达 为 

z -agexp( -bx —cy) (12) 
AP: z =y, E .y-20,.4,:a.b.c 均 为 拟 合 系 
JR. e 是 已 知 的 ,本 刁 不 能 成 为 自 变 量 , 可 以 假定 为 
任意 恒定 值 , 则 拟 合 结果 中 o 不 需要 出 现 ,和 系数 a 
合并 成 一 个 具体 的 数值 。 运 用 最 小 二 乘法 拟 合 图 
10 的 数据 ,得 到 的 模型 方程 如 图 11 所 示 。 由 结 


TF o;, 从 而 在 复杂 三 维 应 力 状态 中 剪 胀 角 模型 写 为 
Y = f(g?,o,,a.) 。 实 验 发 现 水 泥土 单 轴 抗 压强 度 
q, 与 水 泥 掺 入 比 具 有 正 相关 性 ,所 以 可 以 用 gq, 代替 
a。, 则 前 胀 角 模 型 为 Ww = f(g?,o, ,gq,) 。 考 虑 到 前 
胀 角 模 型 方程 两 边 量 纲 的 协调 性 ,将 平均 应 力 和 单 
轴 抗 压强 度 无 量 纲 化 ,并 引入 内 摩擦 角 变 量 , 前 胀 角 
模型 写 为 yy = gf(e?,(g,/g,))。 在 -er 关系 曲 
线 中 ,初始 陡 升 段 的 塑性 应 变相 对 微小 ,忽略 初始 段 
剪 胀 角 并 不 影响 试 样 的 变形 规律 。 剪 胀 角 相 对 于 塑 
性 应 变 及 平均 应 力 是 负增长 关系 , 且 平 均 应 力 与 塑 
性 应 变 取 无 穷 时 剪 胀 角 应 趋 于 无 限 小 。 综 合 以 上 考 


可 知 , 不 同 水 泥 掺 入 比 的 水 泥土 剪 胀 角 模型 方程 拟 
合 相 关系 数 都 较 高 ,反映 出 拟 合 方程 具有 较 高 的 拟 
合 精度 。 拟 合 方程 的 图 形 是 一 个 三 维 曲面 ,在 平面 
上 的 投影 呈 弯 刀 状 ,从 图 11 可 看 出 曲面 能 够 较 好 地 
履 盖 数据 点 。 拟 合 曲面 还 反映 出 剪 胀 角 随 等 效 塑 性 
应 变 增 加 呈 非 线性 下 降 的 特性 ,可 见 剪 胀 角 模 型 方 
程 较 准确 地 描述 了 水 泥土 剪 胀 的 变化 情况 。 事 实 

上 , 因 剪 胀 角 模 型 方程 已 经 隐 含 了 水 泥 掺 入 比 对 前 
胀 角 的 影响 ,从 而 可 以 将 不 同 水 泥 掺 入 比 的 试验 数 
据 进 行 整体 拟 合 ,结果 如 图 12 所 示 , 可见 拟 合 效果 
也 较 好 。 
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Fig.12 Fitting results of dilatancy angle 


of overall cemented soils 
3 KAIRE TEC HU Dr E 


3.1 数值 计算 模型 


为 了 验证 本 研究 所 建立 的 剪 胀 角 演 化 模型 的 合 
理性 ,采用 ABAQUS 软件 对 水 泥土 开展 基于 Mohr- 
Coulomb 强度 理论 和 应 变 软化 模型 的 三 轴 压 缩 数 值 
模拟 试验 ,并 将 数值 模拟 结果 与 三 轴 实 验 结果 进行 
比较 。 

本 次 模型 采用 二 维 轴 对 称 平面 模型 ,模型 尺寸 
与 三 轴 试 验 的 试 样 尺寸 一 致 ,采用 四 边 形 二 次 减 缩 
积分 单元 , 共 划 分 5 个 单元 ,旋转 后 三 维 效果 见 图 
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11 水 泥土 三 轴 试 验 y-e! 关系 拟 合 结 


Fitting results of p-e! relationships of cemented soils using triaxial tests 


13。 模 型 材料 参数 选择 水 泥 掺 入 比 为 40% 的 水 泥 
土 材料 参数 ,弹性 模 量 已 取 各 试 样 的 平均 弹性 模 量 ， 
其 值 为 已 =320 MPa, 泊 松 比 w=0.27, 单 轴 抗 压强 度 
q, =4.0 MPa。 水 泥土 发 生 破坏 后 会 产生 应 变 软化 ， 
所 以 本 次 模拟 采用 的 材料 强度 参数 考虑 软化 效应 ， 
本 构 关系 为 线 弹性 -软化 型 模型 。 


p 1 


到 13 水 泥土 数值 模拟 模型 

Fig: 13 -Cemented soils- numerical simulation model 

在 ABAQUS F, zz fii Hj Mohr-Coulomb 材料 模 
型 , 剪 切 塑性 面 的 硬化 或 软化 是 通过 控制 黏 聚 力 c 
的 大 小 来 实现 。 因 此 ,从 简化 计算 角度 考虑 ,不 考虑 
水 泥土 塑性 硬化 过 程 ,将 黏 聚 力 e 设置 为 线性 下 降 
的 软化 行为 ,如 图 14 所 示 , 而 内 摩擦 角 取 为 恒定 值 
e =43"。 水 泥土 的 剪 胀 角 演化 模型 方程 由 实验 数据 
拟 合 得 到 , 拟 合 方程 见 式 (13 )。 同 时 ,建立 恒定 剪 
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上 胀 角 计算 模型 用 于 比较 ,以 反映 采用 恒定 剪 胀 角 模 
型 进行 计算 所 引起 的 误差 大 小 。 
Y 249. STexp( —15. 22e" -1.4580,/g,) (13) 
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£?/% 


图 14 KIERR IERA 
- Fig. 14 Cohesion model of cemented soils 
SFE ABAQUS 中 ,将 剪 胀 角 演化 模型 方程 通过 
Ffan 语言 写 人 场 变量 子 程序 “USDFLD” 中 ,之 后 
算 模 块 中 调用 子 程序 , 即 可 实现 水 泥土 剪 胀 角 
模型 的 设置 。 
le 
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号 数值 模拟 得 到 的 应 力 - 应 变 、 体 积 应 变 - 轴 向 
关于 曲线 和 实验 数据 进行 对 比 ,如 图 15 ~ 16 所 示 。 
从 图 中 可 以 看 出 ,模拟 结果 与 试验 结果 有 着 较 好 的 
一 鼓 性 。 在 弹性 阶段 ,水 泥土 处 于 弹性 压缩 状态 , 围 
越 大 ,水 泥土 弹性 极限 强度 越 大 ,所 以 体积 压缩 量 
越 大 。 当 水 泥土 破坏 后 ,应力 开 始 下 降 ,而 试 样 体 
积 应 变 开始 反 向 增加 ,增加 的 体积 一 部 分 是 由 于 水 
泥 计 软化 后 回 弹 的 体积 , 另 一 部 分 则 是 由 水 泥土 的 
前 煞 行 为 引起 。 数 值 模 拟 的 水 泥土 剪 胀 量 随 着 塑性 
应 变 的 增加 呈 非 线性 增加 的 趋势 ,增加 的 梯度 在 初 
始 剪 胀 阶段 较 大 ,后 期 逐渐 变 小 。 模 拟 得 到 的 体积 
应 变 变 化 趋势 符合 试验 结果 ,直观 反映 了 所 建立 的 
剪 胀 角 模 型 方程 在 描述 水 泥土 变形 行为 上 的 合理 
性 。 围 奈 是 影响 水 泥土 剪 胀 的 另 一 个 重要 因素 ,在 
低 围 压 下 ,水 泥土 表现 出 较 大 的 体积 膨胀 ,但 随 着 转 
压 的 增加 ,水 泥土 的 剪 胀 行为 变 得 缓慢 。 高 围 压 对 
应 的 体积 膨胀 量 明显 小 于 低 围 压 下 的 体积 膨胀 量 ， 
即 剪 胀 行为 对 低 围 压条 件 更 为 敏感 。 从 模拟 结果 可 
知 ,所 建立 的 剪 胀 角 模型 方程 引入 平均 应 力 的 影响 ， 
对 水 泥土 受 围 压 影响 的 剪 胀 行为 有 着 真实 的 反映 。 

在 以 往 的 工程 计算 分 析 中 ,水 泥土 的 剪 胀 角 通 
常设 为 常数 。 为 了 反映 恒定 剪 胀 角 所 体现 的 剪 胀 行 
为 ,本 研究 设置 了 3 种 恒定 的 水 泥土 前 胀 角 进行 模 


拟 , 分 别 为 yw = eu = pg/2,y = 1°, 对 于 每 一 种 前 
胀 角 取 值 ,在 围 压 设 置 上 选取 了 100 kPa 和 500 kPa 
两 种 工 况 ,结果 如 图 17 所 示 。 从 图 中 可 知 ,在 弹性 
压缩 阶段 ,没有 剪 胀 效应 ,模拟 的 水 泥土 体积 变化 与 
实验 结果 相符 。 当 水 泥土 产生 塑性 变形 后 , 取 值 为 
内 摩擦 角 的 恒定 剪 胀 角 将 引起 试 样 体积 的 急速 膨 
上 胀 , 呈 线性 增长 趋势 ,模拟 得 到 的 体积 应 变 - 轴 向 应 
变 曲 线 严重 偏离 实验 结果 。 而 取 值 为 1" 的 恒定 剪 胀 
角 只 引起 了 微量 的 体积 膨胀 ,并 且 不 足以 抵消 体积 
的 弹性 压缩 ,直至 加 载 结束 试 样 的 体积 变形 仍 是 压 
缩 变 形 ,因此 忽略 剪 胀 角 或 剪 胀 角 取 值 过 小 将 与 实 
际 结果 不 符 。 将 剪 胀 角 取 值 为 o/2 时 ,模拟 结果 与 
实验 结果 相 比 ,模拟 曲线 能 近似 反映 出 水 泥土 在 低 
围 压 区 域 的 体积 膨胀 ,但 水 泥土 在 高 围 压 区 域内 的 
体积 膨胀 被 显著 放大 。 可 以 预知 ,将 剪 胀 角 取 为 其 
他 值 ,也 会 出 现 高 估 或 低估 剪 胀 效应 的 情况 ,说 明 恒 
定 剪 胀 角 模型 不 能 正确 描述 水 泥土 峰 后 的 非 线 性 体 
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Fig.15 Comparison of experimental and simulated 


results of stress-strain curves 
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16 ”体积 应 变 - 轴 向 应 变 曲 线 实验 和 模拟 结果 对 比 
Fig.16 Comparison of experimental and simulated results of 


volumetric-axial strain curves 
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图 17 恒定 前 胀 角 模 型 数值 模拟 结果 


Fig.17 Numerical simulation results with 


constant dilatancy angle model 

剪 胀 角 模 型 的 正确 与 否 ,将 影响 到 水 泥土 材料 
在 剪 切 破坏 过 程 中 的 变形 行为 ,最 终 影响 到 结构 的 
REJ .稳定 性 及 变形 等 工程 问题 。 在 工程 应 用 中 ， 
为 疗 便 计算 常 将 水 泥土 剪 胀 角 取 为 恒定 值 ,这 种 简 
俺 加 作 在 物理 意义 上 并 不 准确 ,之 所 以 采取 这 种 处 
理 方式 ,很 大 程度 上 是 因为 不 清楚 水 泥土 在 峰值 后 
体积 如 何 改 变 52] 。 本 研究 所 建立 的 剪 胀 角 模 型 , 考 
户 平 均 应 力 和 塑性 应 变 因 素 的 影响 ,能 够 合理 表 
述 承 泥土 在 峰值 后 的 非 线性 体积 变化 行为 , 剪 胀 角 
模 济 对 于 指导 水 泥土 在 软 土地 基 加 固 上 的 应 用 具有 
-6 的 工程 意义 。 


> ， 
koi iè 


———— 
的 一 胀 行为 ,建立 一 个 同时 考虑 抗 压强 度 .平均 应 力 
和 驶 性 应 变 影响 的 水 泥土 前 胀 角 模型 ,并 利用 有 限 
元 计算 技术 在 连续 介质 范畴 探讨 了 水 泥土 在 变形 破 
坏 过 程 中 的 体积 变化 行为 ,得 出 结论 如 下 。 

1) 水 泥土 试 样 三 轴 试 验 的 应 力 -应 变 曲 线 呈 软 
化 型 ,在 水 泥土 发 生 破 坏 后 ,伴随 着 明显 的 体积 脱 
胀 。 用 剪 胀 角 来 量化 水 泥土 的 剪 胀 行为 ,结果 表明 
水 泥 斤 入 比 、 围 压 以 及 等 效 塑 性 应 变 会 对 水 泥土 的 
剪 胀 角 产生 较 大 影响 ,实验 现象 反映 出 剪 胀 角 与 水 
泥 挫 入 比 ,应 力 状 态 以 及 等 效 塑 性 应 变 之 间 存 在 着 
客观 上 的 联系 。 

2) 将 围 压 对 水 泥土 剪 胀 效应 的 影响 扩展 到 更 一 
般 的 平均 应 力 的 影响 ,根据 单 轴 抗 压强 度 与 水 泥 摊 
入 比 之 间 的 正 相关 关系 ,用 单 轴 抗 压强 度 代替 水 泥 
挫 和 人 比 的 影响 ,结合 水 泥土 剪 胀 效应 对 等 效 塑性 应 
变 的 依赖 关系 ,建立 了 三 参量 描述 的 剪 胀 角 演 化 模 
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型 方程 。 模 型 方程 拟 合 相关 系数 较 高 , 较 准 确 地 描 
述 了 水 泥土 剪 胀 的 变化 情况 。 

3) 在 ABAQUS 软件 中 对 所 建立 的 剪 胀 角 模型 
开展 三 轴 压 缩 数 值 模拟 验证 ,模拟 结果 证 明 ,考虑 平 
均 应 力 和 等 效 塑性 应 变 影 响 的 剪 胀 角 模 型 能 够 合理 
表达 水 泥土 在 峰 后 的 非 线 性 体积 变化 行为 ,模拟 结 
果 与 试验 结果 有 着 较 好 的 一 致 性 。 而 将 剪 胀 角 取 为 
恒定 值 会 出 现 高 估 或 低估 剪 胀 效应 的 情况 ,不 能 
确 反 映 水 泥土 峰 后 的 非 线 性 体积 变化 规律 。 
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